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ABSTRACT 

Available  world-data  from  flume  experiments  related 
to  bed-load  transport  are  analysed  using  non-dimensional 
Regime  Theory  parameters.  The  analysis  is  based  on  a 
functional  relation  derived  from  the  full' dynamical  state¬ 
ment  of  a  standard  problem  of  bed-load  transport  in  flumes. 

A  relation  between  V^/(s-l)gd  vh’ich  is  a  special  form  of  non- 
dimensional  bed-factor  taking  care  of  buoyant  specific  gravity 

of  bed-load  and  the  relative  deoth  d/  raD  is  proposed  in  the 

v 

d/  mD  range  of  10  to  100  and  charge,  C  range  of  0.0  to  10.0 
v 

parts  per  hundred  thousand.  The  results  obtained  from  the 

analysis  of  flume-data  are  checked  with  a  limited  amount  of 

field  data  and  reasonable  consistency  is  found. 

The  influence  of  other  relevant  non-dimensional 

/  m 

parameters  in  this  range  of  d/  tD  and  C  are  studied.  The 

»  v 

analysis  confirms  the  speculation  that  there  are  definite 
phase  differences  in  different  ranges  of  d/D  and  different 
relations  between  relevant  variables  are  to  be  used  in 
these  different  phases. 

An  annotated  table  showing  the  ranges  of  values  of 
relevant  standardized  parameters,  mostly  non-dimensional, 
for  every  set  of  experiments  is  provided.  Annotation  to 
this  table  is  provided  in  a  parallel  table  shoving  size 
distribution  of  materials  used  in  each  set  of  experiments. 
Brief  descriptions  of  the  shapes  of  these  materials  are  given 
in  various  tables  in  the  main  body  of  the  thesis. 
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CHAPTER  I 


INTRODUCTION 

i.l  GENERAL 

Channels  -  both  natural  and  artificial- 
play  an  important  role  in  the  making  and  unmaking  of 
human  civilization.  Therefore,  man  started,  early, 
to  manipulate  and  construct  them  to  his  advantage. 

But  nature  o'beys  its  own  lavs.  As  a  result,  while 
many  interferences  with  nature  have  been  profitable, 
some  have  caused  wastage  of  money  and  energy  or  even 
disaster,  because  of  lack  of  understanding  of  the  lavs 
of  nature. 

Most  of  the  channels  that  concern  man,  are 
alluvial,  that  is,  they  flow  through  the  materials 
deposited  by  running  water  so  that  their  boundaries 
are  really  self-formed.  Usually  they  still  carry 
the  materials  of  vhich  their  boundaries  are  formed. 

The  science  of  these  alluvial  channels  is 
comparatively  nev  and  a  complete  understanding  of  all 
the  aspects  of  their  behaviour  is  yet  to  be  achieved. 
From  the  end  of  the  last  century,  irrigation  engineers 
in  Indo-Fak  subcontinent  have  tackled  the  problems  of 
alluvial  channels  from  the  viewpoint  of  measured  self- 
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. 
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adjustments  of  artificial  channels  that  refused  to 
retain  their  designed  dimensions.  Relations  were 
developed  for  finally  self -ad justed  channels  with 
marked  success  but  they  could  not  contain  the  measur¬ 
able  factors  of  sediment  load  till  fairly  recently, 
and  then  only  roughly  (Ref.  l). 

In  1909,  Gilbert  (Ref.  2)  took  the  problem 
to  the  laboratory  and  investigated,  to  the  small  scale 
available  there  but  in  terms  of  measurable  bed-load, 
phases  of  transport  outside  of,  and  including  those, 
observed  in  the  field  in  the  former  Indian  continent. 

His  qualitative  observations,  published  in  1914,  added 
to  knowledge  but  his  measured  relations  were  insufficient 
to  produce  an  acceptable  quantitative  theory  for  so 
many  different  phases.  Laboratory  work  has  continued 
till  today  with  widely  contradictory  quantitative  results. 
FIGURES  1.1  to  1.4  from  reference  9  compare  predictions, 
based  on  several  well-known  bed-load  formulas  developed 
mainly  from  laboratory  flume  experiments,  with  the 
"observed"  (measured)  values  in  three  (FIGURES  1.2,  1.3, 
and  1.4)  different  natural  streams.  "Observed"  in 
FIGURE  1.1  means  the  results  obtained  from  modified 
Einstein’s  formula  by  inserting  measured  values  of 
slopes,  velocities,  depths  and  total  bed-material  load 
based  on  suspended  load  samples.  It  Is  evident  from 
these  figures  that  there  is  a  vide  scatter  of  results 
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scatter  as  much  as  SO  to  100  percent  from  their  mean. 

The  authors  (Ref.  3)  admit  that  the  statement  -  "In 
applying  formulas,  errors  of  100  percent  or  more  are 
to  be  expected"  -  by  the  Fontana  Sedimentation 
Conference  1954  (Quoted  from  Ref.  3)  applies  to  these 
cases.  FICrURS  1.5  from  reference  4  adds  a  regime  type 
formula  for  bed-load  transport  in  a  flume,  based  on 
flume  data. 

Till  recently,  mobile  boundary  hydraulics 
has  been  limited  to  the  specialists.  However,  UNESCO 
(NS/NR/17  of  15th  October,  1962)  (Ref.  1)  has  now 
listed  "evolution  of  river-bed  and  sedimentation"  as 
the  third  of  nine  major  scientific  problems  of  hydrology 
for  the  I.H.D.  This  nev  status,  it  is  hoped,  will 
result  in  more  attention  to  the  problems  of  both 
self-adjustment  and  transport  cf  sediment. 

1 • 2  AIM  AND  SCOPE  OF  STUDY 

After  analyzing  the  classic  Gilbert  data, 
Blench  and  Erb  (Ref.  5)  drew  attention  to  the  gaps 
in  knowledge  of  sediment  transport  based  on  flume 
experiments.  They  also  stated  that  d/D  (depth  of  flow 
to  grain  size  ratio)  must  be  relevant  at  certain  stages 
of  flow  and  advocated  the  reanalysis  of  flume  data 
to  examine  the  relevance  of  d/D.  The  present  work  is 
in  line  with  this  thought  and,  here,  an  attempt  will 
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FIGURE  1-5  SEDIMENT  RATING  CURVES  ACCORDING  TO 
SEVERAL  FORMULAS  INCLUDING  A  REGIME  THEORY- ONE 

rased  on  flume  experiments  *  (reprinted  from  ref  * i ) 
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be  made  to  find  this  aspect  of  river-bed  evolution  and 
bed-load  transport.  Therefore ,  the  main  objects  of 
this  study  are: 

1 .  -  to  analyze  available  world-data, 
nearly  all  from  laboratory  flume-experiments 
relevant  to  the  transport  of  bed-materials. 

2 .  -  to  attempt  to  reconcile  the  very 
conflicting  formulas  to-date  and 

3 .  -  to  consider  the  relevance  of  the 
result  to  river-behavio.ur ,  the  aspects  of  which  are 
on  a  far  larger  scale.  The  scope  of  this  study  is, 
however,  confined  to  the  condition  where  bed-load 
transport  is  considered  to  be  very  small. 

Dimensional  analysis,  which  has  been  neg¬ 
lected  till  fairly  recently  (Ref.  5,6,7)  is  used  with 
a  proper  statement  of  the  problem  to  ensure  attention 
to  the  proper  number  of  interrelated  variables.  Regim 
theory  has  been  the  guide  in  choosing  non-dimensional 
parameters  since  it  has  proved  successful  in  the  field 
phase  to  which  it  purports  to  apply  (Ref.  8). 

The  nature  and.  validity  of  the  data  have 
required  study,  through  the  literature,  and  attention 
has  been  given  to  data  on  "initiation  of  motion"  which 
although  different  from  "cessation  of  motion"  -  which 
is  the  bottom  limit  of  transport  -  is  closely  linked 
with  it. 
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1  •  ^  DEFINITIONS  OF  TERMS 

Definitions  of  certain  expressions  and  terms 
used  in  this  study  follow  herein.  Standard  terms  and 
their  symbols  with  dimensions  are  given  in  article  1.4. 
Sediment ;  -  Any  undissolved  material  that  is  transported 
by  a  flow  under  discussion.  Normally,  material, 
denser  than  water  and  carried  by  the  flow  of  water  is 
considered  as  sediment  in  River  hydraulics. 

Bed-load :  -  That  part  of  the  total  sediment  carried 
by  the  flow,  which  moves  by  saltating,  sliding  or 
rolling  along  the  channel  bed. 

Sediment  discharge:  -  The  quantity,  measured  in  units 
of  mass,  that  passes  a  channel  section  per  unit  time 
on  an  average  over  a  period  of  time. 

Bed-material :  -  That  part  of  the  total  sediment  which 
remains  in  touch  with  the  bed  for  most  of  the  times, 
but  may,  at  higher  stages  of  flow,  go  into  suspension. 
Water- sediment  complex : -  The  water  together  with  the 
sediment . 

Charge :  -  Ratio  of  dry  weight  of  bed-load  carried  by 
the  channel  per  unit  time  to  the  weight  of  water-flow 
per  unit  time,  expressed  in  parts  per  hundred  thousand. 
Phase ;  -  State  of  flow. 

Size  or  Grade ;  -  The  size  of  a  sediment  particle  is 
measured  as  the  diameter  of  a  sphere  of  the  same 
material  having  the  same  terminal  velocity  of  settle¬ 


ment  in  water. 


1 


* 

■ 

j 


9 


Critical  Tractive  Force:  -  Tractive  force  (tfdS.)  at 
which  the  general  movement  of  the  bed  starts. 

General  movement:  -  A  condition  of  bed-load  movement 
in  which  bed-particles  up  to  and  including  the  largest 
size  are  in  motion. 

Ripples ;  -  Small  triangular  shaped  bed-forms  that 
appear  on  sand-bed  at  the  beginning  of  motion  of  bed¬ 
load  . 

Bars:  -  Bars  are  deposits  of  bed-materials,  in  the 
channel  section,  that  are  too  low  to  be  classed  as 
islands. 

Dunes:  -  Bed -forms  smaller  than  bars  but  larger  than 
ripples.  In  longitudinal  profile  they  are  triangular 
having  fairly  gentle  upstream  slopes.  They  move 
downstream  and  are  the  principal  sources  of  bed- 
roughness  in  sand-bed  canals. 

Antidunes :  -  Dunes  that  move  upstream  at  supercritical 

flow . 

Regime  theory:  -  Inductive  science  of  channel- 
evolutions  based  on  observations  of  the  dimensions  of 
self -adjusted  canal  systems  in  the  Indo-Gange tic  plains 
(Ref.  1). 

Silt-factor :  -  Lacey's  silt-factor  in  his  basic  regime 
equation,  related  to  silt  size. 

Bed -fa ctcr  and  Side-factor:  -  Lacey's  silt-factor  has 

been  split  by  Blench  into  bed-factor  and  side-factor 

to  account  for  different  natures  of  non- cohesive  bed 
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and  relatively  cohesive  hanks  and  are  defined,  as  V2/d 
and  V^/b  respectively. 

Numeric :  -  A  numeric  is  a  dimensionless  variable. 
Incipient  or  Initiation  condition:  -  Beginning  of  bed¬ 
load  motion. 

Zero  bed-factor:  -  Bed-factor  at  small  bed-load  charge. 
I  A  SYMBOLS  USED  IN  THIS  THESIS 


The  fol loving  symbols  have  been  adopted,  vith 
mass,  length  and  time  as  basic  dimensions.  Other 
symbols  are  described  where  they  appear  first. 

Symbol  Meaning  Dimensions 

A  Cross-sectional  area  of  flow  L2 

a '  Long  axis  of  an  irregularly  shaped 


stone 

b  Breadth  of  the  channel 

b*  Intermediate  axis  of  an  irregularly 

shaped  stone 

C  Bed-load  charge  per  10 

c’  short  axis  of  an  irregularly  shaped 

stone 

d  Depth  of  flow 

D  representative  size  of  bed-material 

J?D  Median  size  of  bed-load  by 

weight 
Bed-factor 

P'bo  Bed-factor  for  small  charge 


L 


L 


L 

L 

L 


L 


L/T 


2 


L/T 


2 
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Symbol 

f 

G 

g 

GR 


\ 


G' 


M 

n 

Q 

q 

R 


S 

s 


Meaning 

Lacey's  silt-factor 

Rate  of  sediment  transport  in 

lbs. /ft.  width 

Body  force  due  to  gravity 

A  representation  of  grain-size 

distribution  which  eouals 


1  (DS4  d50  , 

2  V  '  h6  ] 


Submerged  weight  of  bed-load 
transport  per  sec,  per  foot  width 
Equivalent  sand-grain  roughness 
Roughness  coefficient  of  the  type 
of  Manning ' s  n 

Grain-roughness  coefficient  in 
Meyer-Peter ' s  equation 
Kramer's  uniformity  modulus 
Manning's  roughness  coefficient 
Fluid  discharge-rate 
Fluid-discharge  rate  per  unit 
width  of  channel 
Hydraulic  radius 
Bed-hydraulic  radius  equals 
A/bed  width 
Slope 

Specific  gravity  of  bed-material 


Dimensions 

M/LT 

l/t2 

m/lt 

L 


1?  /t 
l2/t 

L 

L 


_ 

Symbol 

Meaning 
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Dimensions 

t 

Temperature 

© 

T 

Time 

T 

V 

Mean  Velocity  of  flow 

l/t 

vr 

Fall  velocity  of  bed-material. 

l/t 

V* 

Shear  velocity 

l/t 

X 

Non-dimensional  factors  relating 

to  bed-load  size  distribution, 

shapes,  texture  etc. 

gs 

Body  force  per  unit  mass  resolved 

in  the  general  direction  of  flow 

l/t2 

'To 

Tractive  force  at  the  bed 

m/lt2 

To 

Critical  value  of^J  beyond 

which  general  movement  of  bed¬ 
load  occurs 

m/lt2 

IT 

A  constant  -  3.1^ 

- 

y,y- 

Unit  weight  of  fluid 

m/l2t2 

K 

Submerged  unit  weight  of  solid 

m/l2t2 

Unit  weight  of  solid 

.  2  2 
m/l  t 

Mass  density  of  fluid 

m/l5 

fs 

Mass  density  of  bed-material 

M/L5 

1 r 

Kinematic  viscosity  of  fluid 

l2/t 

standard  deviation  in  Charge,  C 

- 

.6 
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The  following  abbreviations  are  used  in  the 
body  of  this  paper. 


UNESCO 

United  Nations  Educational,  Scientific 

and  Cultural  Organisation. 

U.S. 

United  States 

U.S.V.E.S. 

United  States  Waterways  Experiment 

Station 

C.S.U. 

Colorado  State  University 

U.B.C. 

University  of  British  Columbia,  Canada 

Ref. 

Reference 

Fig.,  Figs. 

Figure,  Figures  respectively 

ft. 

Feet  or  Foot  as  the  case  may  be 

R.C.A. 

Research  Council  of  Alberta 

N.R.C. 

National  Research  Council  of  Canada 

F,  -Condition 
bo 

Vanishingly  small  bed-load  charge  - 

Condition 

Indo-Pak 

India  and  Pakistan 

I.H.D. 

International  Hydrologic  Decade 

Indo-Gangetic 

Plains  deposited  by  the  Indus  and 

Plains 

the  Ganges. 

Vs. 

Versus 

P.P.H.T. 

Parts  per  hundred  thousand 

% 

Percent 

B.C. 

British  Columbia,  Canada 

Eq.,  Eqs. 

Equation  and  Equations  respectively 

. 
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gm  . 
cm. 
sec . 

/ 

c  .f . s . ,  cusec. 
mm  . 

fn.  or  fns. 
oC 


Grams 

Centimeter 

Second 

per,  such  as  ft. /sec. 
cubic  feet  per  second, 
millimeter 

function  or  functions 


varies  vith 


CHAPTER  II 


LITERATURE  REVIEW 

2.1  GENERAL 

Before  discussing  the  method  of  analysis  used 
in  this  study,  it  is  necessary  to  discuss  the  sources 
of  data  and  formulas.  The  present  chapter  outlines 
literature  which  is  usually  considered  significant. 

As  experimenters  do  not,  normally,  summarize  the 
ranges  of  important  variables  in  their  data,  and  this 
failure  may  be  relevant  to  the  discrepancies  among 
the  results  (FIGURE  1.1  to  1.5)  the  author  starts  by 
presenting  these  variables.  FIGURES  C-l  to  C-3B 
(Append ix"^  show  the  frequency  distribution  of  sizes 
of  materials  used  by  different  experimenters  whose  data 
have  been  used  in  this  analysis  while  table  B-2 
(Appendix  "B^  lists  the  ranges  of  various  interrelated 
variables  in  each  experimental  set-up  within  the  scope 
of  this  study.  Table  B-2  also  shows  the  ranges  of 
the  values  of  various  standardized  groups,  mostly 
non-dimensional,  for  research  convenience  (Ref.  4). 

2 . 2  REVIEW  OF  LITERATURE  RELATED  TO  MEASURED  TRANSPORT 

DuBoys  (1879)  (Ref.  9)  was  one  of  the  first 
investigators  in  the  field  of  sediment  transport. 
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He  analysed  the  problem  of  transport  by  using  the 
principle  of  conservation  of  energy  and  proposed,  the 
theory  of  "Tractive  force".  His  idea  of  critical 
tractive  force  -  vhich  is  the  tractive  force  at  which 
general  movement  of  bed-materials  takes  place  -  is  still 

shared  by  many  investigators.  DuBoys’  tractive  force 

\ 

equation,  based  on  still  earlier  works  of  DuBuat 
(Ref.  10)  is: 

G  =  Cl^o  (%-%)  (2.2.1) 

where 

G  =  rate  of  sediment  transport  in  lbs. /ft. 
width 

c'l  z  a  coefficient  but  not  a  constant 
^  -  Tractive  force  or  drag  force  at  the 

bed  in  lbs. /ft. 2 
-  Critical  value  of  ^ 

In  1935,  Straub  (Ref.  11)  analysed  the  results  of  various 
investigators  on  the  basis  of  Sq.  2.2.1  and  presented 
the  values  of  c^  and  ^  for  various  sizes  of  sediment 
having  specific  gravity  of  2.65.  These  figures  are 
shown  in  table  2.1. 

TABLE  2 . 1  Values  of  C]_  Calculated  by  Straub 


size  of  sediment 
in  mm.. 

value  of 
cl 

value  of  % 
in  lbs . /f t . ^ 

0.125 

b .  8l 

0.016 

0.250 

0.48 

0.017 

0.50 

0.29 

0.022 

1.00 

0.17 

0.032 

2.00 

0.10 

0.051 

3.0 

0.06 

0.090 
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The  next  and  most  comprehensive  work  in  this 
field  vas  that  of  G.  K.  Gilbert  (Ref.  2)  who 
investigated  the  problem  in  terms  of  measurable 
bed-load,  described  the  various  phases  of  transport 
and  studied  the  individual  effects  of  charge,  flume- 

width,  discharge  and  sediment  size  on  equilibrium  , 

\ 

depth  and  slope.  Gilbert  disagreed  vith  the  commonly 
known  lav  of  sixth  power  (quantity  of  sediment  trans¬ 
ported  varies  vith  the  sixth  power  of  flow  velocity) 
and  mentioned  the  work  of  Deacon  (189^)  who  believed 
that  the  quantity  transported  varies  vith  the  fifth 
power  of  velocity.  Gilbert  explained  that  the  sixth 
power  lav  pertains  to  the  maximum  size  of  the  grain 
that  a  given  current  is  competent  to  move.  This  has 
been  supported  by  William  V.  Rubey  (Ref.  12). 

Engels,  H.  and  Kramer,  H.  (1932)  (Ref.  13) 
presented  a  rather  different  form  of  critical  tractive 
force  equation  in  terms  of  the  physical  characteristics 
of  the  sand  particles  in  a  mixture  -  which  when  con¬ 
verted  to  English  units  reduces  to: 

%  «  0.00138  |  (JT.-Jf)  (2.2.2) 

where 

M  -  Uniformity  modulus  describing  grain  size  . 
distribution 

D  2  Mean  grain  size  in  inches 

<!fi  2  Unit  weight  of  sand  in  lbs. /ft. ^ 

^  =  Unit  weight  of  water  in  lbs./ft.^ 


1  o 
lo 

'T  ~  Critical  Tractive  force  in  lbs. /ft. “ 

Kramer  was  one  of  the  pioneers  to  introduce  size 
distribution  of  the  sediment  as  a  variable  in  the 
problem.  He  described  the  distribution  of  the  sed¬ 
iment  by  vhat  he  called  "uniformity  modulus". 

MacDougall  (1933)  (Ref.  1 4 )  shared  the 

\ 

opinion  about  the  effect  of  grain  size  distribution 
and  pointed  out  the  quite  different  behaviour  of 
uniform  and  graded  sands  though  their  mean  grain  size 
may  be  the  same.  He  agreed  to  the  fact  of  different 
phases  of  bed -movement  but  concluded  from  his  flume 
experiments  that  the  amount  of  sands  moved  is  inde¬ 
pendent  of  the  phases  of  movement  and  that  the  bed- 
load  is  strictly  a  function  of  tractive  force.  He 
proposed  the  following  empirical  relation  for  bed-load 
transport 

G  =  a1Sbl  (Sq-Sqc)  (2.2.3) 

where  G  -  bed -load  in  lbs. /ft.  v id th/sec. 

S  r  slope,  q  -  unit  ^ater  discharge  at  critical 

cond.it  ion 

q  -  water  discharge  in  cusec./ft.  width, 
a-,  are  constants  depending  upon  the  specific 

gravity  and  mechanical  composition  of  the  sand. 

Analysing  mostly  Gilbert  data  supplemented 
by  some  of  his  own,  Schoklitsch  (193*0  put  forward  the 
following  relation  for  uniform  quartz  grains. 

G  -  -3j '•  J-  S  1,5  (Q-bqe) 

D  2 


(2.2.4) 


■ 
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where  o  * 

rate  of  bed-load  movement 

D  ** 

grain  size  in  inches 

S  “ 

slope 

Q  r- 

water  discharge  in  c.f.s. 

b  «  breadth  in  ft. 

q  =  critical  discharge  per  unit  width  of 
\  c 

flume  in  c.f.s. 

From  1932  to  1955  the  U.S.V.E.S.  performed  a 
series  of  experiments  to  determine  the  force  of  flowing 
water  required  to  move  materials  composing  the  bed  of 
the  Lower  Mississippi  River  (Ref.  10).  The  analysis 
of  data  was  based  on  the  DuBoys'  principle  of  tractive 
force  and  the  following  empirical  transport  equation 
was  developed. 


G  = 

1  (dS-doSo)m  (2.2.5) 

n  K1 

where  G  = 

rate  of  movement  in  lbs. /ft.  width/hour 

(dry  weight) 

n  = 

Manning's  roughness  coefficient 

dS  = 

product  of  depth  in  feet  and  slope 

doSo  = 

product  of  depth  in  ft.  and  slope 

at  which  the  linear  plot  of  G  against 

ndSn  meets  the  "dS"  axis. 

m,  are  the 

dimensional  constants  obtained  from  a 

logarithmic  plot  of  the  data.  One  of  the  most  widely 


quoted  criteria  for  the  initiation  of  motion  of 
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non-cohesive  grains  on  flat  bed  is  the  Shields’  diagram 
(FIGURE  2.2.1).  Shields  (1936)  shoved  that  the  tractive 
force  could  be  expressed  as  a  function  of  the  dimen¬ 
sionless  parameter  V* . D/qr  where  D  is  the  diameter  of 
the  sediment  grain.  In  his  experiments  Shields  used 

uniformized  materials  ranging  from  approximately  0.1 

\ 

to  5  mm.  in  size  and  approximately  1 . 06  to  4.25  in  dry 
specific  gravity.  He  concluded  that  non-dimensional 
tractive  force  V  -  now  called  mobility  number  is 
independent  of  V*.  dAt-  now  called  grain  size  Reynolds’ 
number  when  the  latter  exceeds  approximately  500  and 
for  incipient  motion,  with  gravels  exceeding  approx¬ 
imately  3  mm,  in  grain-size,  the  expression  ^ has 
a  constant  value  ranging  from  0.04  to  0,0o.  However, 
his  text  shows  that  he  appreciated  the  possibility  of 
complexities  associated  with  bed-forms.  These  figures 
(0.04-0. Ob ,500)  have  been  disputed  by  Komura  (1963)  and 
Zeller  (1963)  and  Neill  (i960)  (Ref.  15).  Shields, 
apart  from  his  well-known  incipient- condition  relation¬ 
ship,  also  suggested  a  bed-load  formula  (Ref.  16) 


Gs 

'oS 


10 


vhere 


Jf{s-1)D 
q  =  discharge  intensity 
S  =  bed-slope 


(2.2.6) 


s  =  specific  gravity  of  the  sediment 
and  other  terms  are  as  defined  before.  Again  his  text 
shows  he  was  aware  of  possible  qualifications. 


** 
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Shields  diagram  showing  original  data  points 
and  curve  as  usually  quoted. 


Shields-type  curve  by  Iwagaki,  with  grain-size 
values  for  standard-density  material 
as  added  by  Komura. 


Shields-type  diagram  by  Zeller,  with  supplementary  data 
by  Meyer-Peter  and  others. 


FIGURE  2-2-1  SHIELDS'  DIAGRAM  OH  "INI TAT ION  OF 
AND  SHIELDS'  TYPE  DIAGRAM  3Y  IWAGAKI  AND  ZELLER 


MOTION" 
(REF.  15) 


M  -  '  } 
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Vanoni  (Ref.  9)  commented  that  this  formula  applies 
only  to  sediments  which  do  not  form  ripples  or  high 
waves  i.e.  for  sediments  greater  than  approximately  2  mm, 
in  size.  C.  M.  White  (Ref.  17)  derived  two  relations 
for  critical  tractive  force.  For  inviscid  flow  - 
(V,  D/-v  >3.5) 

=  5  (£'*V) D5  tan  q?  (2.2.7) 

O 

where  p  is  a  packing  factor  and  is  equal  to  D  times 
the  number  of  grains  per  unit  area  and  cp  is  the  angle 
of  repose  of  the  material  and  other  terms  are  as  defined 
before.  For  viscid  flow  ( V^D/nr^ 3 . 5 ) 

% 

rX-n-TD  =  C1  p  "5  tanT  (2.2.8) 

6 

where  c^  is  an  experimental  coefficient. 

Experimental  verification  of  Eq.  2.2.8  for  a  50  fold  range 
of  V*D /pr  (from  0.04  to  2.1)  shoved  that  the  value  of  "c^p" 
and  hence  the  value  of  rc/(6-y  w)D  is ’'^dependent  of  V^D/ir 
which  thus  contradicts  the  results  of  Shields  in  this  zone. 
Chang  (1939)  (Ref.  18)  following  and  extending  MacDougall's 
work  proposed  the  following  relation  for  critical  tractive 
force 


where 

sg  =.  effective  specific  gravity  of  submerged 
sand  particles 
c^  -  a  constant 


. 
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D  =  mean  diameter  of  sand  particle  in  mm. 


c 


=  critical  unit  tractive  force  in  lb. /ft. 2 


Using  Sq.  2.2.9  Chang  proposed  a  transport  formula 


of  DuBoys ’  type  for  uniform  sand. 


(2.2.10) 


where  K  is  a  constant  and  all  other  terms  are  as  defined 
before.  He  pointed  out  that  the  advance  of  ripples,  which 
he  assumed  to  be  the  same  as  the  velocity  of  sand  movement, 
is  proportional  to  the  fifth  power  of  the  mean  velocity  of 
the  flow  and  that  irrespective  of  the  modes  of  transport 
sorting  takes  place  in  graded  sands.  One  widely  used  bed¬ 
load  formula  -  that  of  Meyer-Peter  and  Muller  (19^8)  is; 


(2.2.11) 


where  the  term  in  the  square  bracket  is  a  correction  factor 
which  allows  for  wall  effects  and  bed-form  roughness.  In 
Equation  2.2.11, 

G*  =  submerged  weight  of  bed-load  transported  per 


sec.  per  foot  width. 


o  -  ^dS.  where  d  are  unit  weight  and  depth  of 

water  respectively  and  S  is  the  channel  slope. 


-•  mass  density  of  fluid. 


submerged  unit  weight  of  sediment 


D  =  size  of  the  bed-materials 


. 
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q  =.  unit  fluid  discharge,  ^  =  Vmean  x  R^ 

R^r  bed  hydraulic  radius  =  A/bed  width 

^inean  =  ^-b^b  ^ 

=  Grain  roughness  coefficient  =  , 

K  in  feet . 
s 

Equation  2.2.11  is  dimensionally  homogeneous.,.  One 
remarkable  exception  to  this  conventional  approach  of 
tractive  force  theory  to  the  sediment  transport  problem 
is  that  of  H.  A.  Einstein  (19^2)  (Ref.  19).  Einstein 
argued  that  it  was  difficult  to  distinguish,  exactly, 
a  condition  when  bed-materials  start  moving  and  he 
refused  to  use  the  term  critical  tractive  force  for  the 
flow  condition  when  transportation  begins.  He  viewed 
the  transport  problem  as  a  random  process  and  assumed 
that  the  number  of  bed-load  particles  moving  out  of  a 
bed-area  of  unit  width  and  length  per  second  is  a 
function  of  the  probability  that  the  lifting  force  of 
the  flow  exceeds  the  submerged  weight  of  the  particle. 
He  observed  that  the  bed-material  moves  in  steps  or 
jumps  and  assumed  that  the  average  length  of  the  step 
of  an  individual  particle  is  some  function  of  its 
diameter.  Based  on  these  assumptions  and  physical 
reasoning  he  developed  a  relation,  using  dimensional 
analysis,  of  the  form 
G 


m  ) 


(2.2.12) 


(si  -  s)  ■  (  i  >  ai 
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where  G  =  rate  of  sediment  transport/ft.  width 
Vf  s  particle  fall  velocity 
D  ^  diameter  of  grain 


Putting  G/V^D  ,  Eq.  2.2.12  takes  the  form 

4“  fa  (  U-  )  (2.2.12a) 

Einstein,  then  plotted  the  experimental  data  using  Eq. 
2.2.12a.  His  plot  is  shown  in  FIGURE  2.2.2.  Henderson 
(1966)  (Ref.  16)  plotted  the<j> -—relation  in  a  form  due 
to  Brown  (Ref.  20)  who  proposed  the  equation 

4  =  '+0  (  sp-  )3  (2.2.13) 

for  the  upper  straight  line  portion  of  the  curve  shown 
in  FIGURE  2,2.3.  Brown  pointed  out  that  for  low  values 
of  and  hence  of  G,  the  curve  swings  away  from  the 
straight  line  to  the  asymtote  ~  -  0.056  which 
corresponds  to  Sheilds'  initiation  criterion.  Kalinske 
(Ref.  21)  derived  a  relationship  very  similar  to  that 
of  Einstein  hut  he  used  shear  velocity,  V*  in  place 
of  particle  fall  velocity  used  by  Einstein,  Kalinske’ s 
plot  of  G/V* . D  against  Vr  is  shown  in  FIGURE  2.2.4. 

The  constant  value  of  %/r'B  . D  for  initiation  of  motion, 
obtained  from  Kalinske 1 3  plot  is  almost  twice  that  obtained 
from  both  Shields’  and  Einstein’s  plots.  E.  M. 

Laursen  (1953)  (Ref.  22.)  presented  the  following 
relation  for  total  sediment  transport 


. 
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FIGURE  2-2-2(a)  &  (b)  EINSTEIN’S  BED-LOAD  FUNCTION, 
GRAPHS  OF  <J>-y RELATION.  (REPRINTED  FROM  REFERENCE  9). 
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y  (  s  “  1  )D 


FIGURE  2-2-3  HENDERSON’S  PLOT  OF  EINSTEIN'S  BED¬ 
LOAD  FUNCTION  (REPRINTED  FROM  REFE¬ 
RENCE  1 6 ) 


2S 


+  y(S-l)  3) 

FIGURE  2-2-4  KALIN SKE '  S  PLOT  OF  BED-LOAD 


FUNCTION (REPRINTED  FROM  REF.16) 
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C 


(2.2.191  ) 


where  C  =  mean  concentration  of  suspended 


or  bed  or  total  load  as  the  case  may  be 
in  percent  by  weight 


£p  =  sum  of  the  contributions  of  each  fraction 


p  of  size  D. 

shear  stress  due  to  grain  resistance  alone 


-  other  parameters  are  as  defined  before. 

Laursen’s  plot  is  reproduced  in  FIGURE  2.2.5*  Assuming 
that  saltation  accounts  for  all  the  sediment  movement 
Yalin  (Ref.  23)  deduced  a  semi-empirical  relation  for  the 
total  sediment  load.  He  assumed  the  uplift  force  on 
the  particle  to  be  a  function  of  shear  velocity,  V* 
and  from  the  relation  of  uplift  force  and  particle  weight, 
he  worked  out  the  path  of  the  saltating  particle  and 
finally  arrived  at  the  relation 


(2.2.15) 


where  r  -  (%  -%)/%  and 

-0  4 

a  =  2.45.s  (X  /  IT.s.D) 

c 


1 

2 


Yalin’ s  expression  fits  some  experimental  data  for 

/ 

the  particle  size  range  of  0.7  mm.  to  28.6  mm. 
Henderson  (Ref.  16)  described  some  of  Yalin' s  assump¬ 
tions  as  unrealistic  since  the  latter  assumed  that  the 


uplift  force  is  not  dependent  on  time  and  the  increase 


. 


. 

. 
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FIGURE  2-2-5  SEDIMENT-LOAD  FUNCTION  AFTER  E.M. 


LAURSEN ( REPRINTED  FROM  REF . 1 6 ) 
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in  sediment  discharge  is  attributed  to  the  increase 
in  the  range  of  saltation  path  rather  than  to  the 
increase  in  the  number  of  grains  in  motion. 


A 


CHAPTER  III 


THEORETICAL  BACKGROUND  OF  ANALYSIS 


3 . 1  GENERAL 

In  this  chapter  a  brief  history  of  the 
development  and  the  outline  of  the  regime  theory  will 
be  presented  as  the  method  of  analysis  of  this  study 
is  based  on  the  principles  of  this' theory.  A  brief 
outline  of  dimensional  analysis  will  also  be  presented 
here.  Finally,  brief  description  of  the  data  used  in 
this  analysis  will  be  provided.  Detailed  descriptions 
of  these  data,  especially  of  most  experimental  set¬ 
ups  will  be  found  in  reference  28. 

3 . 2  ^REGIME  THEORY  -  HISTORY  OF  DEVELOPMENT  AND  OUTLINE 

The  term  "regime"  has  been  defined  by  Prof. 
Blench  (Ref,  l)  as  "the  behaviour  of  a  channel,  over 
a  period,  based  on  conditions  of  water  and  sediment 
discharge,  breadth,  depth,  slope,  meander  form  and 

progress,  bar  movement,  etc.  - "  and  regime  theory 

studies  the  self-adjustment  of  these  channels  to  long¬ 
term  steady  regime.  Regime  theory  originated  mainly 
from  the  field  experiences  gained  in  construction, 
maintenance  and  extension  of  the  world’s  largest 
canal  systems  in  the  Indo-Pak  sub-continent  since 
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the  latter  part  of  the  nineteenth  century.  The  canal 
systems  in  the  Indo-Gangetic  plains  have  three  degrees 
of  freedom  and  a  tendency  to  acquire  the  fourth  one 
(for  meandering)  and  idealize  the  natural  rivers  so 
veil  that  Prof.  Blench  calls  them  the  "laboratory  of 
the  inductive  science  of  regime  theory/' 

The  idea  of  the  study  of  channel-self- 
ad  .iustment  was  first  conceived  by  R.  G.  Kennedy  in 
1895.  He  studied  the  behaviour  of  a  number  of  self- 
adjusted  canals  in  the  Indo-Gangetic  plains  and  tried 
to  find  a  relation  between  their  hydraulic  and  geometric 
features  and  proposed  the  following  relation 

V  =  0.84d0*64  (3.2.1) 

This  equation  was  one  step  towards  the  design  of  mobile¬ 
boundary  canals  which  used  to  be  designed  by  the 
Iconventional  Manning’s  formula  based  on  rigid-boundary 
hydraulics  but  it  could  not  explain  the  third  degree 
of  freedom  of  the  canal  system  . 

It  was  Lindley,  in  1919,  vho  first  proposed 
the  basic  principles  of  regime  theory  in  its  present 
acceptable  form.  His  statement  was  "when  an  artificial 
channel  is  used  to  convey  silty  water,  both  bed  and 
banks  scour  or  fill,  changing  depth,  gradient,  and 
width,  until  a  state  of  balance  is  attained  at  which 
the  channel  is  said  to  be  in-regime"  (quoted  from 
Blench,  Ref.  l).  Lindley  related  both  bed-width  and 
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depth  vrith  velocity  by  the  following  relations 

V  =  0.95(d)0'5''  (3.2.2.  ) 

V  =  0.57(b)0'535  (3.2.3) 

Gerald  Lacey  (1929,  1933)  proposed  a  change  in  earlier 
regime  relations  and  introduced  his  silt -factor.  He 
proposed  the  use  of  hydraulic  radius  and  vetted 
perimeters  instead  of  Lindley’s  breadth  and  depth 

but  later  Lacey  (1958)  agreed  to  use  breadth  and  depth 
in  sectional  equations.  '  He  also  found  out  a  slope- 
relation  using  his  silt-factor.  Lacey’s  three  basic 
equations  may  be  written: 


Vf 

gtf 


oCf 


1 


R  WsTi/3 

1/6  Q  -  1  A 

vr  '  yp  , _ 

)  X  (p)  xVgRS 


(3.2.4) 

(3.2.5) 

(3.2.6) 


corresponding  to  the  three  degrees  of  self-adjustment . 
N.  K.  Bose  (1936)  proposed  the  use  of  bed-load  size 
instead  of  silt-factor.  Blench  (1941)  proposed 
splitting  of  Lacey’s  silt-factor  to  account  for  the 
different  natures  of  undune d  and  relatively  cohesive 
sides  and  the  non-cohesive  duned  beds  of  mobile¬ 


boundary  channels  and 
side-factor 

2 

f,  =  vr 

b  (i- 


.ntroduced  his  bed-factor  and 


(3.2.7) 


Fs  =  V 


3 


5~ 


(3.2.8) 
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3.3 


He  favoured  the  use  of  breadth  and  depth  in  the 
sectional  equations.  King  (1943)  found,  from  data 
analysis,  a  dynamically  satisfying  slope  equation 

A 


— . —  or  i  Vb  \  ^ 
gdS^  (^=  ) 


(3.2.9) 

Inglis  (19^9)  emphasized  the  inclusion  of  bed-lead 
charge  explicitly  in  the  silt-factor.  Analysing 
classical  Gilbert  data.  Blench  and  Erb  (Ref.  5) 
proposed  a  modified  bed-factor  equation  to  include 


charge  and  thus  converted  King’s  slope  equation 
(Equation  3.2.9)  into  a  sediment  transport  one.  Blench 
and  Qureshi  (1964)  (Ref.  27)  proposed  a  design  curve 
for  rough  estimate  of  in  terms  of  size  of  the 

bed-material,  RJ* and  g. 

BED-LOAD  TRANSPORT  FORMULAS  BASED  ON  REGIME  THEORY 

All  the  bed-load  formulas  and  their  develop¬ 
ment  as  described  in  section  2.2,  with  the  exceptions 
of  Einstein's  and  Kalinske’s  relations  are  based  on  the 
conventional  tractive  force  theory.  A  major  deviation 
from  this  approach  to  the  problem  is  that  of  regime 
theory.  Regime  theory  describing,  as  it  does,  the 
characteristics  of  self -ad justed  channels,  can  be  and 
is  applied  to  the  flumes  which  have  adjusted  themselves 
to  equilibrium  and  the  theory  is,  therefore,  expected 
to  describe  the  transport  of  sediment  as  observed  in 
flume  experiments.  Inglis  (Ref.  24)  made  the  first 


- 
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attempt  to  investigate  the  effect  of  quantity  or 
concentration  of  transported  materials  on  the  regime 
of  the  channels  and  produced  a  relation  to  include 
charge  in  a  regime  equation  as 

bgl/5/Q2  =  fns(Q/D5/2Ei)ni  (CVf/(-*g)1/?)n2 

(5-3.1) 

vhere  C  =  bed-load  charge  in  P.P.H.T. 
b  =  breadth  of  channel 
x  -  kinematic  viscosity  of  fluid 
and  other  terms  are  as  defined  earlier.  Erb  (Ref.  25) 
disputed  the  validity  of  Equation  3*3.1  on  the  ground 
that  it  indicates  the  channels  of  infinite  depth  and 
zero  slope  for  the  important  case  of  vanishingly  small 
charge.  On  the  basis  of  flume  experiments  Blench 
(Ref.  26)  suggested  that  bed-load  charge  be  included 
in  a  bed-factor  equation 

Fb  =  Fbo  U+axC)  (3*3*2) 

vhere  F^  =  bed-factor  =  V~/d 

b 

F  =  bed-factor  for  vanishingly  small  charge 
bo 

and  is  a  function  of  grain  size 
C  =  charge  in  P.P.H.T. 

-  a  constant  . 

After  analysing  Gilbert  and  others’'  data  Blench  and 
Erb  (Ref.  5)  proposed  a  value  of  0.12  for  "a^"  for 
the  condition  of  sub- critical  flow  and  natural  river¬ 
bed  sand.  For  super-critical  flov  a  formula 


■ 
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(3.3.3) 


Fb  -  32.2  +  0.06  (C-Cc) 
was  proposed  for  Gilbert's  uniformised  sand.  Cc  in 
Eq.  3.3.3  is  the  charge  corresponding  to  critical 
velocity.  Blench  and  Erb  also  modified  King's  regime 
slope  equation  into  a  transport  one 

V2/gdS  =3.63  (1  +  C/233)  (Vb/arf  (3.3.4) 
for  sub-critical  flow  only.  Blench  and  Qureshi  (196^) 
(Ref,  27)  suggested,  on  speculative  grounds,  a  relation 

'  Fbo  =  7-3  '  VD  i  *  (■2T0/^)1/6  (3.3.5) 

where  D  -  material  size  in  feet 

-  kinematic  viscosity  at  70°F. 

70 

This  is  believed  to  be  valid  for  all  materials  of 
specific  gravity  2.65  and  larger  than  2  mm.  in  size. 

On  the  same  grounds,  they  also  suggested  a  relation 

V  •  0.5S  V  “**  (*-70«U/re  (3-5.6) 

for  intermediate  sand  range.  The  validity  of 
Equation  3*3.5  is  questioned  by  Neill  (Ref.  15  PP.  27) 

It  is  important  to  note  that  workers  in  regime  theory 
were  aware  that  their  channels  worked  in  a  particular 
phase,'  so  their  formulas  could  not  be  expected  to 
apply  outside  of  some  limits.  However,  non-dimensional 
regime  theory  parameters  may  be  useful  outside  of 
these  limits. 

DIMENSIONAL  ANALYSIS 

"Dimensional  analysis"  .involves  making  a  dimensional 


physical  statement  non-dimensional  to  ensure  the 
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dimensional  homogeneity  that  is  essential  in  the  stat¬ 
ement  of  a  physical  lav.  To  be  of  any  value,  the 
dimensional  statement,  to  which  it  is  applied,  must  be 
the  full  statement  of  a  phenomenon  in  the  form  -  "when 
the  following  quantities  are  specified,  then  the 
following  must  depend  on  them".  The  greatest  advantage 
of  dimensional  analysis  is  that  it  compels  the  full 
physical  statement  of  a  problem.  Its  practical 
advantage-  is  that  it  reduces  the  number  of  independent 
variables  in  a.  statement  to  a  minimum. 

A  'typical  problem  of  alluvial  channels  may 
be  dynamically  stated  in  the  following  words : 

"Suppose  we  have  a  long  straight  horizontal 
flume  of  breadth  b  with  smooth  vertical  sides  at  a 
place  where  body  force  due  to  gravity  is  g.  A  steady 
uniform  discharge  Q,  of  fluid  with  kinematic  viscosity  if 
and  mass  density  flows  through  the  flume  along  with 

sediment,  believed  to  move  as  bed-load,  of  some 

/  m  \ 
representative  size  D  (considered  vD  in  this  paper), 

density  and.  non-dimensional  oroperties  (size 
s 

distribution,  shapes,  their  distribution  and  pattern 
etc.)  X  ,  in  quantity  which  results  in  a  bed-load 
charge  C.  It  might  be  possible  that  the  fluid- 
sediment  complex  is  also  charged  with  suspended 
load  of  specific  size,  shape,  density  and  quantity. 

But  for  simplicity,  these  might  be  excluded  at  this 


- 
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stage.  Nov,  under  these  conditions,  the  channel 
will  adjust  to  equilibrium  vith 


To  the  left  of  Equation  3.4.1,  may  be  added  such 
other  dependent  variables  as  dune-height,  their  wave¬ 
length  and  shear  stress.  But  these  are  not  of 
immediate  interest  in  this  paper.  Dimensional 
analysis  can  reduce  Equation  3.4.1  to 


(3.4.2) 


The  numerics  in  Equation  3.4.2  can  be  commented  as 
follows : 


V  /(s-l)gd  -  is  Blench’ s  bed-factor  amended 


by  (s-l)g  to  allow  for  "buoyant 
weight"  of  sediment.  In  density 


current  works  it  is  known  as 


" dens ime trie  Eroude  Number". 


is  effectively  a  conventional 


d/D 


friction  factor, 
is  depth/grain  size,  so  is  the 
reciprocal  of  what  Neill  (1Qo6) 
calls  relative  grain-size  and  may 


also  be  called  relative  depth  of 


flow . 


G 


is  charge  in  parts  per  hundred 


thousand « 


■ 


4o 


b/d 

f  /f 
's  /  ‘f 

v^/EE-T)  §■ 


is  breadth  to  depth  ratio, 
is  the  specific  gravity  of  bed- 
material  and.  can  be  replaced  by  s. 
is  effectively  a  Reynolds*  Number 
that  Blench  calls  the  VIG  Number . 
Regime  theory  shows  that  it  is 
relevant  in  the  region  where 
(d/D)  >  600  (Ref.  4). 
represents  non-dimensional 
properties  of  the  bed-materials 
as  mentioned  in  the  statement 
of  the  typical  problem. 


3' . 5  SOURCES  AND  DESCRIPTIONS  OF  DATA  USED  IN  THIS  ANALYSIS 

3.5.1  U.S.W.E.S.  : 

These  data  from  the  experiments  carried  out 
by  U.S.W.E.S.  in  Vicksburg,  Mississippi,  were  reported 
in  reference  10.  Nine  different  alleged  natural 
(as  reported  in  the  reference)  sand-mixtures  were 
used  in  the  flume  studies  and  the  rates  of  movement 
of  various  sizes  and  the  tractive  forces  necessary 
for  the  commencement  of  their  movement  were  studied. 

It  was  noted  in  the  paper  that  no  data  were  recorded 
for  any  run  until  an  absolute  equilibrium  had  been 
reached.  But  considering  the  short  duration  of  the 
,  runs  and  considering  the  long  time  which  is  required, 
as  mentioned  by  Blench  (Ref.  2o)  and  Erb  (Ref.  25), 
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by  a  flow  in  flume  to  run  to  regime,  it  was  quite 
probable  that  complete  regime  vss  not  reached  in 
those  experiments.  Besides  the  nine  natural  sand 
mixtures  U.S.W.E.S.  also  used  a  synthetic  sand 
mixture  for  a  series  of  flume-experiments.  The 
synthetic  sand  mixture  was  prepared  by  mixing  together 
uniform  sieve  sizes  of  filter  sand  previously  screened 
from  the  natural  mixtures  of  the  materials.  All  the 
materials  used  in  the  experiments  vere  composed  of 
quartz  and  feldspar  particles  with  small  percentage 
of  other  minerals.  The  mechanical  analysis  curves  of 
these  materials  are  shown  in  FIGURES  C-l  to  C-12.  The 
gradation  Gp  calculated  on  a  method  suggested  by  Simons 
et  al  is  noted  in  tables  3#1  and  3.1(a)  where  other 
physical  characteristics  of  the  materials  are  also 
listed.  The  values  of  D^p,  and 

along  with  the  values  of  Dnn/D  ,  DQ)| /D  and  Dn/D 

^  9o '  50  84'  50  2/  50 

are  shown  in  table  B-l .  Table  B-l  also  lists  the 
values  of  Dgg/D^,  and  D^/D  for  the  Fraser 

River-sand  (B.C.)  for  comparison  with  other  figures 
tabulated  therein. 

In  some  other  experiments  U0S.V.E.S.  (Ref.  37) 
used  light-weight  materials  to  determine  their  suit¬ 
ability  for  use  as  bed-materials  in  the  models. 

Porous  materials  vere  soaked  and  stirred  until  they 
reached  their  maximum  specific  gravity  before 
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Pi  a  cement.  The  general  information  about  these 
materials  is  shown  in  table  3.1(b).  The  grain  size 
distribution  curves  for  each  material  are  shown  in 
FIGURES  C-36  to  C-38. 

3.5.2  T,  Y,  LIU 

T.Y.  Liu  (Ref.  30)  collected  these 
data  at  the  University  of  Iowa  Institute  of  Hydraulic 
Research.  For  incipient  (initiation)  condition, 
depth,  slope  and  rate  of  flow  and  velocity 

/ 

distributions  were  recorded.  The  rate  of  bed-load 
transporation  were  also  measured.  The  material  used, 
was  commercial  sand  collected  from  the  Iowa  River  at 

»  * 

Iowa  city.  This  could  have  been  a  mixture  of  natural 

sands  of  different  D_..  The  grain-size  distribution 

50 

curves  are  shown  in  FIGURES  C-13  and  C~l4  and  the  values 

of  Dc„  D  D„.  ,  D,  and  D0  for  each  material  are 

50'  98'  84*  lo  2 

noted  In  table  B-l.  The  physical  properties  of  the 
materials  are  listed  in  table  3.2. 
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3.5-3  COLORADO  STATE  UNIVERSITY 

Guy,  Simons  and  Richardson  presented  these 
data  in  U.S.G.S.  professional  paper  462-1  (Ref.  31). 

The  experiments  vere  carried  out  in  the  hydraulics 
laboratory  of  Colorado  State  University  during  the 
years  1956-61.  The  purpose  of  the  experiments  was  to 
determine  the  effects  of  the  size  of  the  bed-materials, 

temperature  of  water  and  the  fine  sediment  in  the 

\ 

flow  on  the  transport  phenomenon.  Each  set  of 
experiments  covered  transport  phases  ranging  from  a 
plane  bed  and  no  sediment  movement  to  violent  antidunes. 
The  sands  used  as  bed-material  were  collected  from  five 
different  natural  streams.  The  general  information 
about  the  materials  is  shown  in  Table  3.3  and  the  grain 
size  distribution  curves  are  shown  in  FIGURES  C~15-jt>0 
C-17.  From  the  distribution  curves  shown  for  six 
different  materials,  it  appears  that  only  three  (median 
diameter  0.19,  0.45,  and  0.47  mm.)  materials  follow 
approximately  the  log-normal  distribution  found  in 
natural  river-bed.  The  other  one  (0.27  mm.  median 
dia.)  follows  log-normal  distribution  approximately.  The 

values  of  D^,  D^,  Dp)i,  Dl6  and  D2  for  individual 
material  are  noted  in  Table  3-1.  Photographs  of  each 
sand  showing  their  shapes,  gradation  and  relative  sizes 
are  reproduced  in  FIGURE  3.5.1. 
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FIGURE  3-5-1  PHOTOGRAPHS  SHOWING  SHAPES  OF  SAND  PARTICLES 
USED  IN  C.S.U.  EXPERIMENTS.  THE  NUMBER  IN  THE  LOWER 
RIGHT  CORNER  OF  EACH  PHOTOGRAPH  IS  THE  MEDIAN  DIAMETER 
OF  SAND,  IN  MILLIMETERS. 
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3.5.4  BHARAT  SINGH 

Bharat  Singh  (Ref.  36)  described,  the  studies 
made  (up  to  i960)  in  this  field  by  various  invest¬ 
igators.  In  his  own  flume  experiments  in  the  Hawksley 
hydraulics  laboratory,  London,  he  collected  a  set  of 
data  on  bed-load  transport.  His  flume  was  23  ft.  long 
and  2.5  it.  vide.  Beach  sand  obtained  from  commercial 
suppliers  vas  used  as  bed-material.  The  grain  size 
distribution  curve  from  the  sieve  analysis  of  the 
materials  is  shown  in  FIGURE  C-l8.  Almost  80/  of  the 
sands  vere  between  0.5  to  0.75  nun,  in  size.  The  size 
distribution  curve  is  in  reasonable  agreement  with  the 
log-normal  distribution  found  in  the  natural  river  bed. 
The  specific  gravity  of  the  sand  vas  2.64  and  the 
average  fall  velocity  in  water  at  14°  c.  vas  6.95 
cm. /sec.  The  shape  of  the  grains  vas  fairly  well- 
rounded.  The  time  allowed  to  each  run  to  reach  to 
regime  varied  from  1  to  6  hours.  The  depth  of  flow 
was  adjusted  to  be  uniform  and  the  flow  vas  assumed  to 
be  in  equilibrium.  Regime  conditions  vere  probably 
not  established  yet  the  results  show  a  good  accordance 
with  others  in  this  study. 

3.5.5  WEST  BENGAL  RIVER  RESEARCH  INSTITUTE 

These  data  reported  in  Reference  32  vere 

r 

collected  in  a  tilting  flume  85'  long  and  1.5  ft.  vide. 
The  depth  of  mobile  bed  vas  0.5  ft.  and  the  sides  vere 
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of  glass.  Uniform  flow  was  established  by  adjusting 
the  tail  gate.  Sieved  sand  in  the  range  of  0.15  to 
0.5  aim.  in  diameter  vas  used  as  bed-material.  The 
grain  size  distribution  curve  from  the  mechanical 
analysis  of  the  material  is  shown  in  FIGURE  C-19* 

The  values  of  Dcn,  D  0,  D  ,  D  ,  and  are  noted 

50'  98  84  16  2 

in  Table  B-l. 


5.5.6  BHATTACHARYA 

Bhattacharya  (i960)  (Ref.  55)  carried  out 
flume  experiments  in  the  hydraulics  laboratory  of  the 
University  of  Alberta  to  discover  whether  bed-factor 
,  was  anomalous  when  transport  occurs  in  sub-critical 
flow  without  dunes.  The  experimental  flume  was  12 
inches  wide,  18  feet  long  and  2' -6"  deep.  Experiments 
were  carried  out  both  with  charge  and  without  charge. 
The  range  of  water-discharge  for  runs  without  charge 
was  from  0.52  cusec.  to  1.246  cusecs.  and  the  range 
of  flow  depth  vas  from  2.0  inches  to  11  inches.  Thus 
the  ratio  of  b/d  ranges  from  1.09  to  4.0  only  but, 
of  course,  the  experiments  were  not  expected  to  give 
the  values  of  bed-factor  for  larger  b/d  ratios. 
Considering  these  factors  the  data  were  found  useful 
for  testing  the  effect  of  b/d  on  the  analysis  of  this 
study.  Bhattacharya  considered  most  of  the  runs 

r 

without  charge  as  those  of  F,  -condition  but  some  of 
these  appear  to  be  of  initiation-condition.  The 


. 


. 
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individual  runs  were  allowed  6  to  60  hours  to  have 
the  working  regime  conditions  established.  The 
material  used  as  bed-load  was  1.7  mm.  in  median  diameter 
before  the  experiment  and  1.73  mm.  after  the  experiment. 
The  mechanical  analysis  curve  of  the  material  is  shown 
ip  FIGURE  C-20.  The  values  of  Dcn,  D  D  ,  D  .  and 

5°  98  84  lb 

as  well  as  the  value  of  are  noted  in  Table  B-l. 

3.5.7  DATA  COLLECTED  BY  M.  A.  QURESHI 

Qureshi  (19o2)  (Ref.  34)  virtually  repeated 
Bhattacharya ’ s  experiments  in  a  larger  flume.  The 
material  used  was  the  same  as  that  used  by  Bhattacharya. 

*  The  runs  were  allowed  15  to  72  hours  of  time  to  reach  to 
regime.  The  ranges  of  various  variables  in  Qureshi’ s 
experiments  are  shown  in  Table  B-2. 

3.5.8  OTHER  DATA  ANALYSED  IN  THIS  PAPER 

Flume  data  related  to  bed-load  transport  were 
collected  independently  by  H.  J.  Casey,  M.  P.  O’Brien, 

J.  Bogardi  &  C.  H.  Yen,  Ho  Pang  Yung,  A.  L.  Jorissen, 

C.  H.  MacDougal],  G.  K.  Gilbert,  E.  Meyer-Peter  &  Muller, 
and  Barton  &  Lin.  But  only  small  groups  of  their  data 
fall  within  the  scope  of  this  paper  and  have  been  used 
in  this  analysis.  A  brief  description  of  their  experi¬ 
ments  is  given  in  Table  3.4.  The  mechanical  analysis 
curves  of  the  materials  are  shown  in  FIGURES  C-21  to 
C-33.  The  values  of  D^q,  D9g,  DQip  Dl6,  Dp,  D^/D^ 
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Dq4/d50>  d2/D5o  along  with  the  values  of  Gr  for  each 
of  the  material  used  by  them  are  given  in  Table  B-l. 

A  more  detailed  description  of  their  experimental 
set-ups  will  be  found  in  Reference  25. 

3.5-9  DATA  FROM  INITIATION  OF  MOTION  STUDY  BY  NEILL  (Ref.  40) 

! 

These  experiments  to  determine  velocities  of 
flow  required,  to  cause  "first  displacement"  of  single 
stone  at  various  depths  of  flow  were  carried  out 
in  the  hydraulics  laboratory  of  the  University  of 
Alberta  in  1966.  The  materials  used  in  the  experiments 
were  uniform  in  grain- size  ranging  from  5  to  30  mm.  in 
effective  diameter.  The  flow  condition  vas  uniform  and 
fully  rough.  Considering  the  limited,  width  of  the 
flume  available  the  depth  of  flow  vas  restricted  so 
that  the  b/d  ratio  exceeds  the  roughly  stated  value 
of  5  so  that  there  is  a  mid-central  zone  in  the  flow 
where  the  vertical  velocity  profile  is  the  same  as  in 
an  infinitely  vide  channel.  The  physical  properties 
of  the  materials  are  shown  in  Table  3-5.  FIGURE  3.5.2 
shows  a-  composite  picture  of  various  materials  used 
in  the  experiment.  Neill  also  provided  in  the  same 
report  some  comparable  previous  laboratory  data 
collected,  by  various  investigators.  The  general 
information  about  these  data,  as  provided  by  him. 
is  also  listed  in  the  table. 
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3.5.10  ELBOW-RIVER  DATA 


These  data  reported  by  A.  B.  Hollingshead 
(Ref.  35)  were  collected  in  a  joint  field  project 
sponsored  by  Government  of  Alberta,  Research  Council 
of  Alberta,  University  of  Alberta  and  Water  Survey 
of  Canada.  The  study  reach  was  on  the  Elbow  River  at 
Bragg  Creek  near  C.algary.  Besides  the  general  hydraulic 
survey  of  the  river  including  slope  and  velocity 
observation,  bed-materials  of  the  river  were  sampled 
in  three  different  methods.  The  size  analysis  of 
seven  samples  of  bed-load  are  shown  in  FIGURES  C-3^ 
and  C-35.  The  points  plotted  are  fitted  with  approx¬ 
imate  straight  lines  indicating  a  logarithmic  normal 

distribution.  For  calculation  in  this  analysis  the 
m  "" 

average  vD  was  taken  'as  0.1  ft.  The  values  of  D^g, 

Dg^ ,  and  Dp  for  each  sample  are  noted  in  Table  B-l. 

The  specific  gravity  of  the  material  was  taken  as  2.65 
although  it  was  reported  to  be  2.6k  in  reference  35. 

The  ranges  of  ocher  variables  are  noted  in  Table  B-2. 
The  observed  flow  parameters  are  shown  in  table  D-l. 


CHAPTER  IV 


RESULTS v AND  THEIR  DISCUSSIONS 


4.1  REDUCTION  OF  ANALYSIS  TO  SUIT  DATA  OF  THIS  PAPER 


The  functional  relation 


V2 

'Cs-rjgd 


fns.  (g, 


C, 


8-x))£’£  x) 


fs 


(3.W2) 


shows  that  the  problem  (obviously  a  flume-transport 
one)  is  at  least  seven  dimensional  in  nature  and  the 
analysis  investigating  the  exact  nature  of  relation¬ 
ship  between  the  variables  is  a  complicated  one.  To 
simplify  the  analysis  some  of  the  variables  are  treated 
as  follows: 


1.  C  the  relatively  simple  case  of  low  charge, 
C  is  considered.  Tentatively,  a  value  of  C 
up  to  10  parts  per  1C?  is  considered  as  low 
charge.  To  see  how  C  affects  the  analysis, 
this  range  of  C  (0-10  P.P.H.T.)  is  divided 
into  three  sub-ranges  of  0^0-2  P.P.H.T., 

2/C/5  P.P.H.T.  and  5/C/l0  P.P.H.T. 

2.  b  -  considering  the  general  belief  that  - 
cf 

when  b/d))>5,  the  flow  phenomenon  in  the 
central  zone  represents  that  of  an  infinitely 
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vide  channel  -  the  analysis  in  the  first 

stage  is  confined  to  those  experimental 

data  for  vhich  b /d 5 .  In  the  second  stage, 

to  investigate  the  effect  of  b/d,  all  data 
v 

vith  any  value  of  b/d  are  taken  and  they 
are  divided  into  four  groups  according  to 
their  b/d  values  of  0-2,  2/b/d/5,  5/b/d/l0 
and  b/dj>10 . 

is  essentially  the  specific  gravity  of  the 

solid  and  the  speculation  was  that  /-P 

s7  *f 

might  be  removed  from  its  independent  status 
in  Eouation  3.4.2  if  the  buoyant  weight  of 
bed-load  is  considered  and  g  is  replaced 


by  (s-l)g. 

Equation  3.4.2.  then  reduces  to 


Ti'-rjgd" 


...  fQ 
fns.  (jj. 


'^>r-—~r _ _ _ r  5 

2T  >  fp  X,G^, 


.vhere  Cj_  and.  (b/d)^  nov  represent  certain 
ranges-  of  the  variables  C  and  b/d  respectively. 
From  the  vorld-data  available  upto  da te  mostly  from 
the  collection  of  Prof.  A.  V.  Peterson  and  Mr.  R.  H. 
Cooper  (Ref.  29)  the  parameters  in  Equation  4.1  were 

^  p 

calculated  on  a  IBM-360  computer.  V'/(s-l)gd  is  then 

m. 

plotted  against  d/  t.D  on  a  double-leg  paper  for  the 
data  collected  by  individual  author.  These  plots  of 


_ 
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data  of  individual  author  are  then  superimposed  to 
study  the  overall  nature  of  the  functional  relation 
in  Equation  4.1.  The  plots  are  discussed  in  the 
next  article. 

4 . 2  DISCUSSION  OF  THE  ANALYSIS 

) 

Practical  difficulties  have  prevented  the 

extensive  study  of  all  the  relevant  parameters  in 

Equation  4.1.  But  some  sensible  picture  of  behaviour 

of  V  /(s-l)gd  in  sub-critical  flow  and  within  the 

limitation  of  article  4.1  is  indicated  in  FIGURES  4.1 

to  4.3.  FIGURE  4.1  is  a  plot  of  V2/(s-l)gd  Vs. 

d/  T>  for  a  range  of  C  =  0-2  P.P.H.T.  and  b/d  5. 

It  is  obtained  by  combining  FIGURES  A-l  to  A-o 

(Appendix  “A")  .  FIGURES  A-l  to  A-5  show  the  plots  of 

the  same  parameters  and  for  same  range  of  C  and.  b/d 

as  in  FIGURE  4.1  but,  for  individual  authors'  data 

showing  points  representing  various  materials  used 

in  the  experiments.  FIGURE  A- 6  is  a  plot  of  several 

.  authors'  data  named  in  the  figure  itself.  Although 

the  plotted  points  in  FIGURE  4.1  show  some  scatter, 

they  indicate  a  general  trend  in  the  behaviour  of 
2 

V  /(s-l)gd  within  the  reasonable  width  of  a  band. 

This  trend  is  discussed  in  the  latter  part  of  this 
article . 

FIGURE  4.2  is  a  plot  of  V2/(s-l)gd  Vs. 
d/  D  for  a  range  of  2/C/5  P.P.H.T.  and  b/d  >  5. 
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FIG.  4-3  PLOT  OF  V7(s-1)  gd  vs.  d/ 


. 


This  is  again  obtained  by  combining  FIGURES  A-7 
to  A-ll  vhich  are  in  turn  plots  of  individual  authors’ 
data  for  the  same  range  of  C  and  b/d.  FIGURE  4.7 

o  m 

shows  the  plot  of  V  /(s-l)gd  Vs.  d/  D  for  5/C/lO 
P.P.H.T.  and  b/d  ^5  and  is  obtained  by  combining 

FIGURES  A-12  to  A-lo  which  are  again  counterparts  of 

) 

FIGURES  A-7  to  A-ll  but  with  5/c/lO.  To  study  how 

m  -7  _ _ r\ 

material  sizes  (  p)  and  for  that  matter  V2/g 

P  m  *^r 

affect  the  relation  between  V  /(s-l)gd  and.  d/  D, 

v 

FIGURES  4.4  (a)  and  4.4  (b)  were  plotted  with  data 

m 

subdivided  into  various  ranges  of  D  values  as  noted 

v 

in  these  figures.  FIGURES  4.5  (a)  and  4.5  (b)  are 
plots  of  similar  type  but  with  data  subdivided  into 
different  ranges  of  Gpj  values.  FIGURE  4.5  (c)  is  a 
plot  of  H.  J.  Casey  data  who  used  two  materials  of 
quite  different  gradation  G^. 

Scatter  of  points  in  FIGURES  4.4  (a)  to  4.5  (c) 
prevents  a  definite  statement  about  the  dependence  of 

P 

V~ / ( s  —  1 ) g d  on  the  par  ticle  sizes,  and  their  distribution. 

FIGURE  4.6  shows  bands  containing  about  95/ 
of  the  points  of  FIGURES  4.1,  2,  7.  -The  flattened  bottoms 
of  the  bands  suit  known  behaviour  of  V*/:(s-l  )gd  in  irri¬ 
gation  canals  and  are  consist ant  with  these  points 

o 

although  not  proved.  V  /(s-l)gd  apparently  depends  on 
C,  perhaps  in  proportion  to  about  (l  0.10C). 
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The  composite  hand  ABCD  in  FIGURE  4.6 

encloses  the  plotted  data  of  FIGURE  4.1  to  4.5  where 

charge  C  varies  from  0.0  to  10.00  P.P.H.T.  This  hand 

m 

has  a  slope  of  -4/5  in  the  d/  ,D  range  of  10  to  100. 

v 

The  central  line  of  this  band  appears  to  follow  the 
relation 


(s'-l)gd 


(4.2.1) 


m  m 

within  this  range  of  d/  D.  Beyond  d/  D  value 

v  w 

greater  than,  approximately  400,  the  lines  limiting 

the  hands  tend  to  become  horizontal  shoving  the 

p  m 

independence  of  V“/(s-l)gd.  on  d/  t  D  in  this  zone. 

From  observations  in  the  Indo-Gangetic  canal  systems 

where  "P  /-P  3nd  X  factors  are  practically  fixed  by 
s  »f 

nature,  regime  theory  workers  explain  that  beyond 

m  p 

d/  0^400  (approximately),  V  /(s-l)gd  is  dependent 
v 

on  l/ifg(s^TJ  D/iT.  This  is  also  noticed  in  FIGURES 

4.4  (a)  and  4.4  (b)  where  points  representing  larger 

size  of  materials  fall  mostly  above  those  representing 

the  smaller  size  materials. 

m 

Between  d/  D  values  of  100  to  400,  there 
is  curious  transition  between  those  two  phases 
described  earlier.  This  transition  seems  to  be 
parallel  to  that  found  in  conventional  friction 
factor  diagram  and  in  this  transition  zone 


o  U1  ■?  _ 

V  /(s-l)gd  appears  to  depend  on  d/  D,  C ,  yi^/ST s ~  1 )  D/lT 

and  m3y  be  on  "X"  factors  too.  Sufficient  flume-data 
m 

covering  d/  D  range  of  1  to  10  were  not  available. 

V 

As  a  result,  the  nature  of  relation  between  the  relevant 
parameters  in  this  zone  could  not  be  reliably  studied. 
Moreover,  when  V  /(s-l)gd  exceeds  about  0.6  the  flow 
becomes  supercritical.  Available  bed-load  transport- 
data  from  flume  experiments  with  low  charge  (C/10 
P.P.H.T.)  and  covering  this  phase  (supercritical) 
of  flow  are  insufficient.  But  from  the  limited  amount 
of  transport  data  plotted  in  FIGURE  4.7  and  from  the 
fair  number  of  "initiation"  data  collected  by  Neill 
and  others  and  plotted  in  FIGURE  4.8  some  speculations 

r\ 

may  be  made.  These  "initiation"  data  cover  a  V  /(s-l)gd 

m 

range  of  0.6  to  2.0  and  d/  D  range  of  2  to  10.  When 

v 

FIGURE  4.6  is  superimposed  on  FIGURES  4.7  and  4.8 

p 

it  appears  that  the  relation  between  V^/(s-l)gd  and 

n  m 

d/  yD  in  the  d/  t  D  range  of  2  to  10  follow  closely 

m 

their  relation  in  the  d/  D  range  of  10  to  100. 

v 

The  three  bands  in  FIGURE  4.6  each  represent¬ 
ing  one  of  the  three  ranges  of  C  as  noted  earlier 
show  to  a  certain  extent  that  as  C  increases  the 

band  as  a  whole  moves  upward  indicating  a  slightly 

2  m 

higher  value  of  V  /(s-l)gd  for  same  d/  D  with  an 
increase  in  C.  This  gives  the  indication  that  the 
scatter  of  plotted  points  within  the  band  itself 
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may  be  due  to  the  variation  in  C  within  the  range  .of 

0-10  P.P.H.T.  To  verify  the  effect  of  C  within  the 

fitting  band  of  FIGURE  4.6,  y~/(s-l)gd  was  plotted 

against  d/  D  in  FIGURE  4.9  where  the  plotted  data 

were  split  into  three  ranges  of  C  that  is  C=0.01-0.2, 

0=1. 8-2.0  and  0=9.0-9-99  P.P.H.T.  It  is  evident 

from  FIGURE  4.9  that  points  representing  higher  values 

of  C  fall  above  those  representing  lower  values  of  C 

indicating  that  the  scatter  within  the  band  is  due  to 

the  variation  in  charge  within  its  range  of  0-10  P.P.H.T. 

2  m 

FIGURE  4.9  (a)  shows  a  plot  V  /(s-l)gd  Vs.  d/  *D  for 

w 

data  on  light  weight  materials  with  various  specific 
gravities  mentioned  in  the  figure.  These  experiments  ver 
done  by  U.S.V.S.S,  to  examine  the  suitability  of  these 

materials  for  use  as  bed-materials  in  models.  In  the 

\ 

same  figure  Neill’s  "initiation"  data  with  cellulose 
acetate  (specific  gravity  -  1.31)  are  also  plotted. 

When  FIGURE  4.9  (a)  is  studied  along  with  FIGURE  4.6  it 
appears  that  the  plotted  points  in  figure  4.9  (a)  fall 
within  the  band  of  Figure  4.6  in  spite  of  their  low 
specific  gravities.  Also,  no  systematic  effect  of 
specific  gravity  can  be  noticed  on  the  scatter  in  FIGURE 
4.9  (a).  This  leads  to  the  belief  that  the  parameter 
f./fr  in  Equation  4.1  can  be  taken  care  of,  if  data  are 
plotted  in  terms  of  V2/(s-l)gd. 
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FIGURES4.10  to  4.1 3  are  the  results  of  an  attempt 

to  investigate  the  effect  of  b/d  on  the  plot  of 
o  m 

V  / ( s-1 )gd  Vs.  d/  D.  FIGURE  4.10  is  a  plot  of 

V  /(s-l)gd  Vs.  d/  D  for  G  =  0-2  P.P.H.T.~and 

v 

FIGURES  4.11  and  4.12  are  plots  of  same  parameters 
as  FIGURE  4.10  but  for  2/C/5  and  5/C/lO  P.P.H.T. 
respectively.  Points  plotted  in  FIGURES  4.10  to 
4.12  are  grouped,  irrespective  of  their  authors, as  -per 
ious  ranges  of  b/d  values  as  mentioned  in  article 
4.1  and  also  noted  in  these  figures.  These  plots 
(FIGURES  4.10,  4.11  and  4.12)  are  obtained  by 
combining  FIGURES  A-18  to  A -21  ,  A-22  to  A-2'5  and 
A-26a  to  A-29  respectively.  FIGURES  A-18  to  A-29 

p  m 

are  plots  of  V  /(s-l)gd  Vs.  d/  D  for  individual 

v 

authors'  data  and  for  the  different  ranges  of  b/d 

and  C  as  noted  earlier. 

FIGURES  4.10  to  4.12  show,  on  the  whole, 

the  same  kind  of  relation  between  V  /(s-l)gd  and 
m  m 

d/  D  in  the  range  of  d/  D  between  10-100  as  found 

in  FIGURES  4.1  to  4.9  and  discussed  earlier.  The 

intermixing  in  the  position  of  the  plotted  points 

irrespective  of  their  four  different  ranges  of  b/d 

indicates  that  b/d  does  not  have  any  systematic  effect 

m 

on  the  plot  within  the  d/  range  of  10  to  about 

500. 

FIGURE  4.13  is  a  counterpart  of  FIGURE  4.6 
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v.ith  b/d  range  of  about  1.JL0  to  120.  When  FIGURE 
4.15  is  compared  with  FIGURE  4.6,  where  the  plotted 
data  come  from  the  experiments  where  b/d  ^5,  no 
difference  can  be  found  between  the  two  and  this 
leads  to  the  belief  that  the  plot  is  independent 
of  b/d  within  the  ranges  mentioned  earlier. 

FIGURES  4.14  and  4.15  are  plots  V2/(s-l)gd 
Vs.  d/  for  the  data  collected  by  Bhattacharya 
(Ref.  55)  and  Qursshi  (Ref.  54)  respectively.  The 
b/d  ratio  in  Bhattacharya 1 s  experiments  were  all 
less  than  5  and  the  range  of  b/d  values  in  Qureshi's 

\ 

experiments  varied  from  1.95  to  7.5.  In  both  the 
cases  the  bed- condition  was  reported  to  be  that 
of  F^0  -  condition  in  most  experimental  runs.  In 
spite  of  their  low  values  of  b/d,  the  points  fall 
fairly  well  within  the  band  of  FIGURE  4.15. 

\ 

4  •  ^  TREATMENT  OF  "INIATION  OF  MOTION"  DATA 

FIGURE  4.16  is  a  plot  of  V2/(s-l)gd  Vs. 

d/  D  for  the  data  collected  by  Neill  (i960)  to 

study  the  first  displacement  of  bed-materials. 

FIGURE  4.17  is  a  similar  plot  of  the  comparable 

data  of  other  investigators  provided  by  Neill. 

FIGURES  4.l6  and  4.17  plotted  together  result  in 

FIGURE  4.18  which  shows  close  accordance  with  the 

m 

results  discussed  in  Article  4.2  in  the  d/  D  range 

w 

If  FIGURE  4.18  is  superimposed  on 


of  10  to  100. 
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FIGURE  4.6,  it  will  be  seen  that  most  of  the  initiation 
data  fall  close  to  the  bottom  limit  of  the  band  for 
C  =  0-2  P.P.H.T.  The  positions  in  FIGURE  4,l6  of 
plotted  points  pertaining  to  light  weight  materials 
(cellulose  acetate)  whose  specific  gravity  is  1.31 
indicate  that  the  parameter  *f  / -f  in  Equation  4.1 
may  be  effectively  taken  care  of  by  considering  the 
buoyant  weight  of  sediment  and  by  replacing  g  by 

(s-l)g. 

4 . 4  TEST  OF  ANALYSIS  WITH  FIELD  DATA 

Field  data  observed  on  the  Elbow  River 

described  in  Article  3-5.10  were  available  to  test  the 

results  obtained  by  analysing  flume-data.  FIGURE  4.19 

shows  a  plot  of  these  data  along  with  data  on  gravel- 

bed  irrigation  canals  in  the  San  Luis  Valley  (Colorado) 

collected  by  Lane  and  Carlson  (Ref.  38). 

Table  D-l  and  D-2  show  their  observed  flow 

m 

parameters.  Although  the  data  cover  a  d/  D  range 

v 

of  about  10  to  60  the  plot,  within  this  range,  looks 

i. 

consistent  with  the  results  of  flume  experiments. 

4  •  5  COMPARISON  WITH  BED-LOAD  DATA  ANALYSIS  BY  J.  ROTTNER 

J.  Rottner  (Ref.  6)  analysed  about  2,500 
flume  data  recorded  by  different  investigators  and 
compiled  by  J.  V.  Johnson  (Ref.  28).  Applying 
dimensional  analysis  to  the  problem  he  came  up 
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with  the  following  non-dimensional  parameters 

V 


a  densimetric  Froude 

number  i 

where  G  is  solid  discharge  per 
unit  width 

as  defined  before 


VPTs^rj 

g _ 

fVgd^(s-l)  — 

d 
D 

q 

—  which  is  again  a  kind  of 
Froude  number  _  G 

Rottner  then  plotted  y=  V/x/oXs-'^T  against  x  -jf^/gd  ( s - 1 ) 
for  various  ranges  of  d/D  values.  In  these  parameters, 

V  is  in  cm. /sec.,  G  is  in  gm./sec./cm.  width,  d  is  in 
cm.  The  exponent  l/3  was  chosen  for  the  graphs  in 
order  to  render  the  various  degrees  of  scatter  com¬ 
parable.  A  sample  of  Rottner’ s  plot  is  shown  in  FIGURE 
A-3  0.  From  these  plots  he  developed  an  equation  of  the 
form 


1/3 


x  =  axy  -  bx 


(4.5.1 ) 


where  x  and  y  are  as  defined  before  and  a^  and  b^  are 
some  constants.  Now  for  charge  C-^0  the  value  of  G  in 
Rottner ’s  parameter  tends  to  zero.  Therefore'  x->0 . 
Therefore  from  Equation  4.5.1,  for  incipient  condition 


y  - 


V 


or 


v/arsmj  = 


‘i 


0.5.2) 
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2 

Squaring  and  dividing  both  sides  by  g  in  cm. /sec.  . 

We  get 

V2  ,bl  2  1 

WW-T7  -  (5J)  g  (4.5.?) 

The  left  hand  side  of  Equation  4.5.3  is  the  parameter 
plotted  against  d/D  in  this  analysis. 

From  his  equation 

x  -  axy  -  bx 

Rottner  calculated  the  optimum  values  for  the  co¬ 
efficients  a^  and  b^  by  using  the  method  of  least 
squares  and  from  this  he  deduced  the  characteristic 
b^/a^  for  incipient  entrainment  condition.  His 
tabulated  values  are  reproduced  in  Table  4-1.  From 

o 

these  tabulated  values,  V  /gd(s-l)  values  were 

calculated  by  using  the  Equation  4.5.3  taking  g 

as  981  cm. /sec. In  FIGURE  4.20  the  plot  of  V“/gd(s-l) 

against  d/D  for'  these  values  is  shown.  The  curve 

is  in  full  agreement  vith  other  plots  of  this  analysis 

in  the  d/D  range  of  10-100.  In  his  paper  Rottner 

2 

remarked  that  the  slope  of  the  plot  V  /(s-l)gd  Vs, 
d/D  is  not  constant  but  seems  to  increase  with  d/D. 

This  is  also  noticed  in  FIGURE  4.20. 
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TABLE  4-1 


TABLE  SHOWING  cl/D 

AND  b1/a1  VALDES 

ADOPTED 

FROM 

ROTTNER  (1959) 

Page 

Mean.  (d/D) 

(b1/a1) 

1 

6.71 

22.5 

2 

8.38 

19.6 

3 

10.1 

19.-4 

4 

12.9 

17.6 

5 

15.3 

16.5 

6 

17.3 

14.6 

7 

20.1 

14.8 

8 

22.6 

13.7 

9 

25.6 

13.7 

10 

28.3 

11.9 

11 

31.3 

10.5 

12 

34.0 

11.6 

13 

36.7 

11.0 

14 

39.7 

9.69 

15 

45.0 

10.4 

16 

46.6 

10.1 

17 

49.8 

8.01 

18 

5^.2 

8.97 

19 

58.1 

8.82 

20 

61.5 

8.49 

21 

65.3 

7.48 

• 
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TABLE  4-1  CONTINUED 


Page 

Mean  (d/D) 

<Vai) 

22 

'  69.0 

8.03 

23 

71.9 

7.52 

24 

75.7 

8.15 

25 

79.6 

6.95 

26 

84.5 

7.39 

27 

89.1 

6.96 

28 

93.7 

7.3 

29 

99.7 

7.37 

30 

105 

7.00 

31 

112 

6.66 

32 

119 

OJ 

% 

VO 

33 

126 

6 . 08 

3^ 

134 

5.77 

35 

142 

5.14 

36 

151 

5.75 

37 

lol 

4.83 

38 

178 

4.92 

39 

193 

4.05 

40 

215 

4.38 

41 

233 

4.26 

42 

262 

3.90 

43 

298 

4.46 

V 


, 


TABLE  4-1  CONTINUED 

Page 

Mean  (d/D) 

(^l/al ) 

44 

550 

5.06 

45 

472 

5.56 

46 

870 

3.51 

* 


NOTE : 


"Page”  in  Rottner's  analysis  stands  for  various 
plots  he  classified  according  to  various  ranges 
of  d/D. 
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CHAPTER  V 


CONCLUSIONS  AND  RECOMMENDATIONS 


5.1 


CONCLUSIONS 

The  following  conclusions  can  be  drawn 
from  the  results  obtained  by  analysing  flume-data 
and  their  reasonable  consistency  with  the  limited 
amount  of  field  observations: 


a) 


(ii) 


(iii) 


(iv) 


In  the  range  of  data  analysed  the  behaviour 
2 

of  V  /(s-l)gd  is  quite  different  from  that 
in  the  range  of  sand-bed  canals  and  rivers. 
The  data  in  the  range  10/d/  ^  D/lOO  fall  in 
a  band  whose  centre-line  can  be  represented 


by: 


linger  =  (d/D)"7 


(5.1) 


The  band  limits  (Line  AB  and  CD  in  FIGURE 
4-6)  containing  about  95 /  of  points  can  be 
represented  reasonably  by  replacing  the 
coefficient  3.7  in  Equation  5.1  by  8.3  and 
1.8  respectively. 

The  regime  theory  relation 


2 

cs-hp 


fn  (l^gTs'-TT^j 


(5.2) 
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r 
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(vi) 


(vii) 


(viii) 


can  be  seen  coming  into  effect  above  values 
m 

of  d/  D  in  the  neighbourhood  of  400. 
v 

The  change  from  Equation  5.1  to  Equation 
5.2  is  fairly  abrupt  and  seems. to  be  of 
different  types;  data  are  insufficient  to 
permit  a  definite  statement  on  the  types 
and  to  reach  appreciably  into  the  zone 
of  Equation  5.2. 

The  literature  does  not  describe  bid-forms 

sufficiently  to  show  a  reason  for  the  phase 

change  from  Equation  5.1  to  Equation  5.2. 

The  scatter  of  points  about  equation  5.1 

does  not  appear  to  be  due  to  the  variations 

m 

in  D;  but  805>  of  the  data  had  D  values 

v 

between  0.2  mm.  and  2.0  mm.  so  a  slow 

variation  of  V^/(s-l)gd  with  D  could  be 

hidden.  But  a  variation  of  V2/(s-l)gd  with 

C  Is  remarkably  noticed. 

A  systematic  effect  of  b/d  could  not  be 

found  for  1.10  /  b/d  /  120  approximately; 

neither  grain-size  distribution  expressed  by 

Gp,,  ranging  from  1.07  to  4.5,  seems  to  have 

m 

any  systematic  effect  in  the  range  10  /  d/  D 

v 

m 

l  100.  But  in  the  range  100  /  d/  D  /500, 

v 

data  with  smaller  values  of  G^  seem  to  fall 
above  those  with  larger  values  of  G^>  (FIGURE 

4.5  (a)  &  4.5  (b)). 


. 
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(ix) 


Light  weight  materials  plotted  consistently 
with  the  others  in  terms  of  V“/(s-l)gd. 


5.2 


RE  C  OMMS  EDAT I Oh S 

The  following  points  merit  attention  for 
future  research: 

(i)  The  relatively  blank  range  1  /  d/D  /  10 
should  be  investigated . in  order  to  make 
definite  statement  about  the  behaviour  of 
V^/(s-l)gd  in  this  zone. • 

(ii)  Experiments  should  be  repeated  with: 

a)  River-bed  sand  of  standardised  grading 

b)  Coarser  materials  manufactured  to  have 
the  standard  river  sand  grading 

c)  d/D  reaching  to  about  1,500  so  as  to 
cover  the  start  of  the  regime  theory 
phase . 

(iii)  Special  experiments  with  artificial  particles 
of  fixed  sizes  and  shapes  (e.g.  spheres  etc.) 
and  with  various  specific  gravities  are 
recommended . 


(iv) 


Data  covering  higher  amount  of  bed -load 


charge  should  be  analysed. 
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